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Dissipative quantum control of a spin chain 
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A protocol is discussed for preparing a spin chain in a generic many-body state in the asymptotic 
limit of tailored non-unitary dynamics. The dynamics reqnire the spectral resolution of the target 
state, optimized coherent pulses, engineered dissipation, and feedback. As an example, we discuss 
the preparation of an entangled antiferromagnetic state, and argue that the procedure can be applied 
to chains of trapped ions or Rydberg atoms. 


The robust generation of entangled states is a corner¬ 
stone in quantum technological applications. Multipar¬ 
tite entanglement, in fact, plays a crucial role in various 
tasks of quantum information and communication 
and can enable one to achieve unprecedented levels of 
precision in sensor systems and metrology 0. High- 
fidelity generation of multi-partite entangled states has 
been demonstrated with protocols based on deterministic 
operations [7]. These protocols realize dynamics which 
optimally couple an initial state to the target state with 
a sequence of unitary operations. Their performance be¬ 
comes more challenging as the number of components N 
increases. In fact, this often implies a larger number of 
high-fidelity operations, which makes the protocol more 
sensitive to parameter fluctuations and to disorder, and 
requires longer time scales, over which the detrimental 
effects of intrinsic noise and decoherence become more 
relevant. This situation has motivated the search for al¬ 
ternative strategies. 

One promising approach is based on engineering noise 
[S] and dissipation [SHU] in order to drive a many-body 
system towards the desired non-classical target state. 
Protocols based on this idea are often denoted by quan¬ 
tum reservoir engineering (QRE), and their hallmark is 
the robustness against parameter fluctuations, which re¬ 
sults from the non-unitary nature of the processes that 
pump the system into the target state. When based on 
dissipation, they can be considered a many-body gen¬ 
eralization of optical pumping, originally proposed by 
Kastler for creating spin polarized atomic ensembles by 
means of spontaneous decay [13] . As in optical pumping, 
the target state is stable, effectively decoupled from the 
mechanism which pumps out all other states involved 
in the dynamics, but fed by the dissipative processes 
[nidsi- Under these premises the population of the tar¬ 
get state will increase asymptotically towards unity. 

The formal procedure for implementing QRE is usually 
based on constructing a Liouvillean £ for the density ma¬ 
trix p of the system, for which the target state qt is the 
unique stationary state, i.e. Lqt = 0 [HI [T6|. When 
Qt = IV't) (V'tI, then the condition can often be cast in 
terms of rate equations, which couple the population of 


the target state, Pp = Pt:{qt p{ty\, with the populations 
Pn of the states |V'n), forming together with \'tpT) a com¬ 
plete and orthogonal basis in the state space. Denoting 
by Fn-^m > 0 the rate coefficients for the transitions 
\ipn) —t IV'm)) the equation for Pp reads 

Py = -FpPr + Y, ^n^rPu , (1) 

n^T 

and the loss rate of state \ipn) is r„ = The 

objective is to achieve Pp —)■ 1 as t —> oo. From con¬ 
siderations based on detailed balance, it can be veri¬ 
fied that a necessary condition for efficient production 
of pt is r„ ^ Er for n ^ T. It is sufficient when 
min„^ 7 ’ „i (Tn-j-m) ^ Tr, which can be reached by ex¬ 
ploiting symmetries of the dynamics m- This idea is 
at the basis of several proposals for dissipatively pump¬ 
ing spin or harmonic-oscillator systems into bi-partite 
and into specific multi-partite entangled states; examples 
are found in Refs. [ll[igilMlj- Experimental demon¬ 
strations include realizations with trapped ions |27I129j . 
atomic ensembles [3D], and superconducting qubits m- 
The identification of the procedure, however, becomes 
more complex for arbitrary multi-partite entangled tar¬ 
get states. This calls for the development of viable proto¬ 
cols for the dissipative preparation of a generic entangled 
state of a many-body system. 

Here we discuss such a procedure. Our idea is to tailor 
the excitation spectrum of the many-body system, such 
that the target state is an eigenstate and all transitions 
between pairs of states can ideally be individually ad¬ 
dressed. Engineered dissipation allows one to perform 
irreversible population transfer, in order to construct dy¬ 
namics as in Eq. ([^. Transitions from the target state 
are far off resonance from all pumping processes, so that 
the outcoupling rate Ft is sufficiently small. By repeat¬ 
edly applying a sequence of pulses that empty all other 
states, the system is pumped into the target state, with 
an asymptotic fidelity that depends on the ability to tai¬ 
lor the transition rates. This procedure generalizes a 
method for quantum-state preparation of molecules [32] 
to many-body systems. It provides a complementary ap¬ 
proach to the one proposed in Ref. |33j . In addition, in 
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order to counteract noise and decoherence, which become 
more and more important as the number of components 
increases, measurements followed by feedback operations 
are built into the pulse sequence which restore the effec¬ 
tiveness of the procedure for long evolution times. 

We illustrate the procedure by discussing the prepara¬ 
tion of a spin chain in the entangled antiferromagnetic 
state. 



(b) 


E 


iV’T) = (inn..-) + unt-.-))/v^, ( 2 ) 

where | f) and | |) are the two energy eigenstates of a 
(pseudo) spin 1/2, separated by huj^. For N = 2 ions, 
is a triplet (Dicke) state, which can be perfectly 
decoupled from collective spin excitations via quantum 
interference processes [34]. One can thus construct dy¬ 
namics for which IV't) (or the corresponding singlet state) 
is stationary; examples are in Refs. millllSillSSj- For 
N > 2, however, this procedure cannot be directly ap¬ 
plied, since |'0 t) is no longer a Dicke state. 

In order to realize non-unitary dynamics of which \^pT) 
is the stationary state, one could construct a harmonic 
Hamiltonian of which {iPt) is the ground state. Dissipa¬ 
tive preparation into the target state would then proceed 
by means of a generalization of sideband cooling |37] . Im¬ 
plementing these dynamics with spins would in general 
require one to work with an equidistant energy spectrum 
whose excitations are collective spin states, thus in gen¬ 
eral Dicke states. It has the drawback that, as N grows, 
the number of undesired ” dark states” that are decoupled 
via quantum interference increases and that disorder and 
inhomogeneities may render the spectrum anharmonic. 

Instead, we choose a Hamiltonian whose spectrum is 
purposely tailored to be anharmonic and of which |'0r) is 
an eigenstate, but not necessarily the ground state. The 
interaction Hamiltonian 

(3) 

i=i 

serves this purpose, with Jj > 0 and the Pauli op¬ 
erator for spin j. This Hamiltonian can be realized for 
chains of Rydberg atoms [3H] or of trapped ions (35] . In 
the latter case it is implemented by tailoring the coupling 
between the internal and the external degrees of freedom 
of the chain: the coupling coefficients Jj depend on the 
interparticle distances and are symmetric about the cen¬ 
ter: in a linear Paul trap Jn-i = Ji and Ji < Ji+i for 
j = 1,..., [^]. Fig. j^a) shows the specific case = 4. 

Hamiltonian (|^ stabilizes the target state and iden¬ 
tifies the states V'n entering the rate equations, whose 
coefficients Tn^m shall be engineered. The transitions 
\tpn) Itpm) are driven resonantly by laser pulses, 
whereby the detunings A from the spin transitions vary 
from pulse to pulse. The corresponding spin Hamiltonian 
in the rotating frame reads 

Ho(A) = -AA^a|/2, (4) 
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FIG. 1: Procedure for preparing the spins of a four-ion chain 
in the entangled state (| t4t4) + I 4t4t))/'\/2 via non-unitary 
dynamics, (a) The spins interact with one another {Jj) and 
with an ancilla, here a harmonic oscillator of frequency ujt . (b) 
The energy spectrum of the spin chain is tailored by control¬ 
ling the couplings Jj (only the In = 0) manifold of the ancilla 
is shown). Arrows indicate the pulses which resonantly pump 
the spins into the target state; each arrow represents two op¬ 
erations, as illustrated in (c): one coherent pulse (solid) that 
entangles spins and ancilla, and one dissipative pulse (dashed) 
that makes the transfer from state to state irreversible. In 
(d) the pulse frequencies corresponding to the arrows in (b) 
are displayed, using the same colour code. A denotes the 
detuning from ujq. The dashed arrows in (b) indicate transi¬ 
tions driven by the ’’blue” and ’’green” pulses that resonantly 
pump the spins into a state that is not the target. This state 
is depopulated by the ’’grey” pulse. 
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of which \'4’t) is eigenstate for N even [301 • The resonant 
transitions for = 4 are shown in Fig. 

The desired asymmetry in the coefficients and 

Fm-^n is achieved by means of engineered dissipation, 
along the lines of Ref. m- It is realised by first en¬ 
tangling a single spin with an ancilla, followed by dis¬ 
sipation on the ancilla (see Fig. [^a,c)). In our ex¬ 
ample the ancilla is a harmonic oscillator, which may 
be a lossy cavity mode coupling with the spins |2Ij . 
or a collective vibrational mode of an ion chain, cou¬ 
pling with the spins via the mechanical effects of light 
and sympathetically cooled to the ground state as in 
Ref. |5H1. Denoting by |0) and |1) the oscillator’s ground 
and first excited state, irreversible population transfer 
from |' 0 „) to IV’m) is performed by first coherently driv¬ 
ing the transition |■^/'n)| 0 ) —t |^/>m)|l), and then damping 
the mode by an external reservoir that induces the tran¬ 
sition |V’m)|l) ^ |V’m)|0). 

The coherent laser-driven dynamics which entangle 
spins and ancilla are described by the Hamiltonian 

H{A) = Hl^^+Ho{A)+h'^gjaj{a+a^)+huJta'^a. ( 5 ) 
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Here, gj is the Rabi frequency, whose value is sufficiently 
small in order to drive only resonant transitions, and a 
is the annihilation operator of the harmonic oscillator 
at frequency w*. The oscillator is cooled at rate 7 to a 
steady-state excitation number fi <C 1 ; the non-unitary 
cooling dynamics are described by the superoperator [IT] 

-\-l)'D[a]p + ^nD[a)]p, ( 6 ) 

where 'D[X] is a functional of the operator X such that 
'D[X]p = XpX^ — {X^Xp + pX^X)I2, and p is the den¬ 
sity matrix of spins and ancilla. Pumping into the target 
state is realized by sequences of pulses, whose compo¬ 
nents correspond to the map T(A,t) = ^ 

which alternates Liouvillian £coh(A)p = [H{A), py{ih) 
for time t with engineered dissipation as in Eq. (IM for 
time y. The protocol iterates the concatenated map 

T = T{Aj, tj)'r{Aj-i, tj-i)... 7”(Ai, ti) , (7) 

where the sequence, the detunings and the durations tj 
are optimized to drive the system into the desired asymp¬ 
totic state with close-to-unit fidelity. Over the time scale 
of a sequence, the dynamics can be cast in terms of a rate 
equation as in Eq. Q and the choice of the detunings re¬ 
sults in tailoring the effective coefficients. The idea might 
be regarded as a dissipative extension of the Law-Eberly 
protocol [ 32 ] , originally developed for arbitrary quantum 
state preparation by means of coherent dynamics, and 
based on identifying the individual steps which determin¬ 
istically connect an initial and a final state. Indeed, with 
our procedure, we achieve pt = lim^_>.oo Tran{T^p( 0 )} 
for a certain set of initial states p( 0 ), where Tran de¬ 
notes the trace over the ancilla’s degrees of freedom. 
Nevertheless, we find that the pumping efficiency con¬ 
siderably drops if the initial state is an equal statistical 


mixture and/or if the Rabi frequencies gj vary signif¬ 
icantly. In fact, in these cases pumping happens into 
both {iPt) and the degenerate antisymmetric superpo¬ 
sition (I titi ■ • ~ I itit ■ ■ ■)) /Vs. One remedy could 

be to alternate Hamiltonian H{A) with another pump¬ 
ing Hamiltonian, assuming one can engineer the spatial 
gradient of the pulse phase, but this approach is not ro¬ 
bust against decoherence and fluctuations in the values 
of the couplings gj. 

Our solution that enables pumping from arbitrary ini¬ 
tial states into the target state is to include in the se¬ 
quence a parity-correcting operation based on the proto¬ 
col of Refs. 1331131] ■ It performs a parity measurement, 
described by the operator 

11 = ay., ay, ( 8 ) 

followed by a conditional operation on the system which 
corrects the parity in the case that it is not the desired 
one. The corresponding dynamics can be realized by 
means of an ancilla, whose relevant states are denoted 
by |H) and \B): given the state of the system is \tp) and 
the ancilla is prepared in |H), first the unitary map 

\y)\A) ^ 

is applied; this map is the identity if the state possesses 
even parity, while if the state has no dehnite parity it 
becomes entangled with the ancilla. If in a subsequent 
measurement the ancilla is found in \B), conditional dy¬ 
namics are performed that invert the parity of the sys¬ 
tem’s state; in our simulations, a az operation is applied 
to one of the spins. Another option could be to reinitial¬ 
ize the spins to | f • • ■ t) • Denoting by V the correspond¬ 
ing map, the complete sequence of pulses we implement 
is y = VT. This protocol is efficient for arbitrary initial 
states and has constant depth, which makes it scalable 
[45] . It is conceptually an extension of methods for cool¬ 
ing the motion of ions based on measurements |3SHlHj 
(see also Ref. [33] for an application to QRE). Its realiza¬ 
tion requires that the coupling between spins and ancilla 
is homogeneous to a good degree. 

We now come to the specific features of an imple¬ 
mentation based on Hamiltonian (|^, for the example 
of the target state of Eq. ([^ with = 2,4 ions in a 
linear Paul trap. We first identify the frequencies, i.e. 
the detunings Aj of the pulses, for all transitions which 
couple the target state to any other state. For 2 ions, 
these are Ai = —cot + 2Ji and A 2 = — Ai, which pump 
I tt)|0) -t |V’t)| 1) and | U)|0) |V’t)| 1), respectively; 

afterwards the ancilla is damped, making the transfer 
irreversible. The Rabi frequencies of the pulses are set 
to spectrally resolve the individual resonances; residual 
off-resonant coupling, which would depopulate the target 
state, is minimized by choosing pulse durations for which 
this coupling produces an integer number of Rabi oscil¬ 
lations. For this choice, we pump N = 2 ions into the 
target state with fidelity X > 1 — I0“® under ideal con¬ 
ditions, i.e., for h = 0 and in absence of other sources of 
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FIG. 2: Residual infidelity for preparing N = 4 spins in the 
target state of Eq. it for n = 0 and in absence of deco¬ 
herence. The sequence T' contains pulses at the detunings 
Ai = —Wt + 2(Ji -I- J 2 ) (red arrows in Fig. |^b)), A 2 = 
-LJt + 2Ji (blue), A 3 = -ujt - 2(Ji - J 2 ) (greyT A 4 = -Ai 
(green), and A 5 = — A 2 (purple). The other parameters are 
Ji = J 3 = 0.05u;t, J 2 = V^Ji, Qj = g = 5 X 7^7 = 20. 

The pulse durations are such that Jitj > 1 {Jitj ~ 10), the 
pulse areas are optimized to minimize the loss rate Ft [49j . 
The procedure is robust against parameter fluctuations: the 
infidelity doubles when the pulse areas change by about 30%. 


noise. For TV = 4 ions we identify a sequence of 5 pulses, 
shown in Fig. [^b) and detailed in the caption of Fig. 
which leads to an asymptotic fidelity > 0.9995 under 
ideal conditions. The infidelity as a function of time is 
displayed in Fig.[^ 

We now analyse how the fidelity is affected by the abil¬ 
ity to engineer the desired dissipation. We first consider 
varying the temperature of the reservoir and thus h. Fig¬ 
ure a) displays the asymptotic fidelity for different val¬ 
ues of n and for = 2,4 ions, and shows that the control 
on engineered dissipation becomes more stringent as the 
number of spins is increased. Figure |^b) shows the fi¬ 
delity in presence of noise and decoherence, which we 
consider here to be due to spin flips at rate yflip and 
dephasing with yueph- The corresponding Liouvillians 

are yaip Y,j (t’K] -b T>[cr]]^ and yueph Y.j and are 

added to £coh(^)- The parity operation V counteracts 
the decoherence and keeps the fidelity above 0.9 for values 
of yflip, ydeph for which in absence of V we observe a drop 
to 0.5. The curves show that the effect of noise becomes 
more detrimental as N grows, which is expected as the 
protocol becomes slower because of the spectral crowd¬ 
ing around the target state. This could be counteracted 
by increasing the energy splittings in the spectrum, here 
given by Wj and Jj. It is important to note that the 
strength of the coupling, determining the speed of each 
pulse, scales differently with N depending on the phys¬ 
ical system. For instance, when the dissipative channel 
is the collective motion of an ion chain, the coupling be¬ 
tween spins and motion decreases as N grows, due to the 




FIG. 3: Fidelity for preparing N = 2 (blue squares) and 
N = 4 (red circles) spins in the target state as a function of 
(a) n (in absence of decoherence of the internal state) and (b) 
the rate of incoherent processes 2 yflip in units of g (for h = 0 
and assuming yuip = yseph). The fidelity corresponds to the 
asymptotic value of a sequence of pulses with the map T' ■ 
The fidelity is independent of the initial state, provided the 
parity correction is included in the sequence. The parameters 
are given in the caption of Fig. 

increasing inertia of the crystal. If instead it is a cav¬ 
ity mode, the coupling can increase with -x/fV owing to 
superradiant emission. 

In conclusion, we have described a procedure for 
preparing a spin chain in an entangled state. The pro- 
cotol uses engineered dissipation, which makes it robust 
against moderate fluctuations in the parameters, and in¬ 
cludes a parity-error correction procedure, which makes 
it robust against detrimental noise and decoherence. The 
basic requirement is that the target state is spectrally re¬ 
solved, which is achieved by constructing a suitable spin- 
spin interaction. Efficient preparation of iV-spin entan¬ 
gled states is warranted as long as the required spectral 
resolution is larger than the typical rate of noise and de¬ 
coherence. For arbitrary initial states, the protocol time 
scale is expected to increase exponentially with N, as it 
requires the capability to sweep over all state space. It 
can be notably reduced if the initial state is, e.g., a po¬ 
larized chain such as | f, t: ■ ■ • i t)i which is typically easy 
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to produce with optical pumping. Its duration can be 
further shortened by optimizing the time duration of the 
coherent pulses, e.g., by using time-dependent values of 
A and g in Eq. (§ identified by means of optimal-control 
theory [50] . 
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